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Database of optical properties of cosmic dust analogues (DOP) 
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Abstract. We present a database containing information on different aspects of calculation and usage of the optical 
properties of small non-spherical particles — cosmic dust analogues. The main parts of the DOP are a review of 
available methods of the light scattering theory, collection of light scattering codes, special computational tools, 
a graphic library of the optical properties (efficiencies, albedo, asymmetry factors, etc.), a database of the optical 
constants for astronomy, a bibliographic database of light scattering works, and links to related Internet resources. 
The general purpose of the DOP having the address http://www.astro.spbu.ru/DOP is twofold — to help scientists 
to apply the light scattering theory in astronomy and to give students and beginners a possibility quickly to get 
necessary knowledge on the subject. 
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1. Introduction 

Dust has been detected in most astronomical objects - 
from the Solar system to active galactic nuclei, and every- 
where the dust grains were found to play an important role 
in different physical and chemical processes. Modelling 
of these objects should include accurate consideration of 
scattering, absorption and emission of radiation by the 
dust grains. 

So far in such a modelling one mainly utilized the 
model of homogeneous spherical particles, though there 
are many evidences that the interstellar grains are mostly 
non-spherical and inhomogeneous. Today enough possibil- 
ities exist to avoid the use of the spherical model. To help 
scientists to make this step, we have developed a database 
including useful information, data, codes, references and 
links and made it accessible via the Internet. As far as we 
know there are no analogous databases in the World Wide 
Web. 

In this paper our database having the address 
http://www.astro.spbu.ru/DOP is described first briefly, 
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and then partly in more detail. The database contains var- 
ious introductory notes written for students and beginners 
and more specific information on light scattering theory 
and its applications. The DOP includes many pieces of the 
review "Optics of cosmic dust" by Voshchinnikov ( 2002 ) 
and materials being used in lectures on physics of inter- 
stellar matter in the St. Petersburg University. 

2. General overview of the database 

The information, resources, and links presented in the 
DOP are collected in three main blocks: 

1) those related with the theoretical aspects — scat- 
terer models, methods of their realization, bibliography; 

2) those connected with calculations (and hence ap- 
plications) — optical constants, light scattering codes, 
graphic and tabular libraries of the optical properties, spe- 
cial tools; 

3) those coupled with other topics (e.g., radiative 
transfer modelling). 

The database has the following sections and subsec- 
tions: 
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A. Theoretical aspects 

Al. Definitions of the optical characteristics of non- 
spherical scatterers 
A2. Scatterer models and their parameters 
A3. Exact and approximate methods of light scattering 
theory 

A4. Light scattering experiments 

A5. Bibliography of light scattering works 

B. Application aspects 

Bl. Optical constants of materials 
B2. Light scattering codes 
B3. Benchmark results 

B4. A graphic library of the optical properties 
B5. Self-training algorithm for calculation of the optical 
properties of fractal-like aggregates 

C. Related resources 

CI. Radiative transfer tools 
C2. Miscellaneous 

Let us now turn to a brief description of the subsections. 

2.1. DOP subsections on theoretical aspects (part A) 

Here we first define the cross-sections, albedo, scatter- 
ing matrix and other basic optical characteristics of non- 
spherical scatterers and compare them with those for 
spherical particles. The idea is to supplement and com- 
press information given in other sources of such definitions 
(see Mishchenko et al., 12000112113) 

The parameters of non-spherical scatterers (size, 
shape, etc.) are briefly discussed in the next subsection 
and the light scattering models popular in astronomy 
(spheres, cylinders, spheroids, aggregates) are listed in 
there. 

To realize any scattering model, i.e. to obtain its op- 
tical properties, one should select a proper exact or ap- 
proximate light scattering method. We make an overview 
of various light scattering methods and cite several other 
useful reviews of this kind. A special very detailed re- 
view of the methods was prepared by N.G. Khlebtsov (see 
Sect. 3.1 for more details) 

The optical properties of cosmic dust analogues can be 
also measured in laboratory. References to two available 
WWW databases of such measurements are presented. 

The last subsection is devoted to bibliography on light 
scattering theory and its various applications which is con- 
sidered in detail in the Sect. 3.2. 

2.2. DOP subsections on application aspects (part B) 

Here we start with definitions of optical constants ( refrac- 
tive indices or dielectric functions) of materials and link 
the Jena-St. Petersburg Database of Optical Constants for 
astronomy (JPDOC) developed earlier. This database de- 
scribed by Henning et al. (1999) and Jager et al. (2003) 
consists of original laboratory data, references and links to 
resources related to measurements and calculations of the 
optical constants for materials important for astronomical 



applications. More information on this part is presented 
in Sect. 3.3. 

The next subsection contains a set of light scattering 
codes developed by us and used in astronomical applica- 
tions. Links to other collections of such codes are also 
given. 

To test a code not used before one often needs 
some well checked results. The optical characteristics of 
spheroidal particles of different size, shape, and refrac- 
tive index calculated by several methods are given in our 
benchmark subsection. 

An important part of the database is the graphic li- 
brary of the optical properties of different scatterer models 
- homogeneous and core-mantle spheres, homogeneous 
and core-mantle cylinders, homogeneous spheroids. In an- 
other part of this subsection the properties of dust grains 
which can be estimated from the observations of the in- 
terstellar extinction and polarization are discussed (see 
Sect. 3.4 for more details). 

The DOP also includes description and computer re- 
alization of a self-training algorithm to predict the optical 
properties of fractal-like aggregates. As this tool was not 
yet properly outlined in the literature, Sect. 3.5 is devoted 
to its brief description. 

2.3. DOP subsections with related resources (part C) 

The optical properties are very often used in radiative 
transfer (RT) calculations. Therefore, we concern this 
point as well as list the RT methods and the majority 
of the RT codes used in astronomy. 

The last subsection contains links to different related 
resources (like other databases of light scattering codes, 
electronic newsletters, etc.) which might be of some use in 
the DOP context. 

3. A more detailed description of some parts of 
the database 

3.1. Review of light scattering methods for 
non-spherical particles (A3) 

The review consists of two parts — an bibliographic 
overview of history and development of several basic meth- 
ods and a comparison of all available approaches to solu- 
tion of the light scattering problem for non-spherical par- 
ticles. 

The overview made by N.G. Khlebtsov concerns both 
the exact (separation of variables, point matching, integral 
equations, coupled dipoles, T-matrix) and approximate 
[Rayleigh, Rayleigh-Gans-Debye (RGD) and its general- 
izations, anomalous diffraction, eikonal, geometrical op- 
tics (GO), perturbation, etc.] methods. In this review the 
history of development of the methods is traced in very 
careful manner. The works on the applicability range and 
various applications of the methods are referred as well, 
and all together there are over 500 cited papers. 
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Fig. 1. Applicability ranges (domains left to the curves) of 
different methods in the case of dielectric prolate spheroids 
(a/b is the semiaxis ratio, xy the size parameter of equiv- 
olume sphere). Calculations were made with codes based 
on the Rayleigh (RA) and quasistatic (QSA) approxima- 
tions, the discrete dipole (DDA), separation of variables 
(SVM) and T-matrix (TMM) methods. 



The comparison of methods is made in form of a table 
for 11 main approaches: 

— separation of variables method, 

— discretized Mie formalism, 

— finite element method, 

— finite difference method, 

— point matching method, 

— generalized multipole technique, 

— extended boundary conditions method, 

— method of moments, 

— coupled dipoles method, 

— Fredholm integral equation method, 

— ray tracing/Monte Carlo method. 

These approaches are classified and their ability to 
treat particles of different shape and structure as well as 
anisotropic and chiral scatterers and ensembles of particles 
is considered. For each method, references to the ansatz 
and some recent reviews are given. Comments to the table 
show the relationships and differences of the approaches, 
briefly concern their classifications and features, sources 
of available codes, and refer to works on comparison of 
the methods and other reviews. 

We have also started the numerical investigation of 
applicability ranges of most popular methods using the 
spheroidal model of scatterers. An example of the results 
is given in Fig. 1 (see Il'in et al., 20(21 for more details). 

3.2. Bibliography on light scattering and its 
applications (A5) 

Besides the review, the DOP also includes a special search- 
able bibliographic database on light scattering theory and 



its various applications. This bibliography being collected 
for over 15 years by V.A. Babcnko consists of the follow- 
ing parts: 

— spherical homogeneous particles — methodological 
questions, functions, calculations, internal fields, reso- 
nances, near field, etc. 

— inhomogeneous spheres — layered, with radial and 
chaotic internal inhomogeneities, with inclusions, reso- 
nances, etc. 

approximations for homogeneous spheres - - GO, 
Rayleigh, RGD, eikonal, van de Hulst (anomalous diffrac- 
tion), and so on approximations; 

— complicated problems — beam scattering, scattering by 
charged spheres, by rough, anisotropic and chiral parti- 
cles, etc.; 

— cylindrical particles (see below) ; 

— non-spherical objects — works on different approaches 
and methods; 

— specific scattering media — interstellar particles, ma- 
rine particles, metallic zoles, powders, etc.; 

— carbon; 

— ice; 

— water; 

— agglomerates, fractals; 

— non-linear and mechanical effects in aerosols; 

— vaporization and explosion of particles; 

— inverse problems; 

— laboratory experiments; 

— laser remote sensing; 

— radiative transfer, 

— etc. 

To demonstrate the wideness of the questions consid- 
ered we shall exemplify a full list of subsections for the 
section on cylindrical particles: 

— circular infinite cylinders; 

— oblique incidence and orientation averaging; 

— multi-layered cylinders; 

— GO approximation (including hexagonal cylinders); 

— cylinders with smooth radial inhomogeneity; 

— gyrotropic, gyroelectric and anisotropic cylinders; 

— internal and near fields; 

— experiments on cylinders; 

— morphological resonances; 

— van de Hulst approximation for cylinders; 

— Rayleigh and RGD approximations; 

— Hart-Montroll and S-approximations; 

— non-circular infinite cylinders; 

— circular cylinder of finite length and discs. 

The total number of included references exceeds 8000 
and each subsection contains several dozens of references. 

3.3. Optical constants for cosmic dust analogues (Bl) 

This part of the DOP includes some notes introducing the 
optical constants required by astronomy, and links the op- 
tical constants database (JPDOC). The JPDOC consists 
of two main parts. One is the bibliography on optical con- 
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Table 1. Jena data currently available from the JPDOC 



Compound 


Composition 


State 


Low temperature 
measurements 


Spectral re 


gions 


silicates 


(Mg,Fe)Si0 3 , (Mg,Fe) 2 Si0 4 


glassy 


yes 


ultraviolet 


- infrared 




MgSiOs, (Mg,Fe) 2 Si0 4 


crystalline 


yes 


infrared 






MgSi^Oy 


amorphous 


yes 


ultraviolet 


- infrared 




(Ca,Al,Mg,Fe)Si a; !/ 


amorphous 


110 


infrared 




sulfides 


(Mg,Fe)S 


crystalline 


yes 


infrared 






SiS 2 


crystalline 


110 


infrared 




oxides 


(Mg,Fe)0 


crystalline 


yes 


ultraviolet 


- infrared 




AI2O3 


amorphous 


no 


infrared 






(Mg,Al)0, 


crystalline 


no 


infrared 




carbon 


a-C:H 


amorphous 


no 


ultraviolet 


- infrared 



stants of astronomically interesting materials: 

■ amorphous/glassy/crystalline silicates of different 
kinds; 

- silicon, SiO, crystalline/fused Si0 2 ; 
— metals: Fe, Mg, and others; 

- oxides: FeO, Fe 2 3 , Fe 3 4 , MgO, A1 2 3 , MgAl 2 4 ; 

- sulfides: FeS, MgS, SiS 2 ; 

- carbides: SiC, FeC, TiC; 

- carbonaceous species: diamonds, graphite, coals, kero- 
gens, HAC, glassy /amorphous carbon, PAHs and so on; 

- organics: tholin, "organic refractory", etc. 

- ices: H 2 0, CO, C0 2 , NH 3 , HCN, etc. and their mix- 
tures; 

- FeSi, CaC0 3 and some other materials. 

We tried to include all papers published since the begin- 
ning of the previous century (all together over 700 refer- 
ences) . 

Another part is a set of data specially measured in 
the Jena laboratory for astronomical purposes. Note that 
although various terrestrial analogues of cosmic solids 
have been studied in chemical and physical laboratories, 
many of these experiments neither took into account the 
specifics of cosmic dust materials (composition, lattice 
structure, etc.) and conditions (low temperature, process- 
ing, etc.), nor covered the wavelength intervals of the cur- 
rent astrophysical interest. 

A summary of the data obtained in Jena and made 
available is given in Tabled data examples are presented 
in Fig. 2 (see Jager et al., 2003 for more examples). 

3.4. Graphic library of optical properties of 
non-spherical scatterers (B4) 

This DOP section consists of two parts — a basic graphic 
library and special notes on optical properties of cosmic 
dust analogues and their connection with observed mani- 
festations of interstellar dust. 

In the basic library we show how different optical char- 
acteristics (extinction, scattering and absorption efficien- 
cies, albedo, asymmetry factor, etc.) vary with scatterer 
parameters: the real (n) and imaginary (k) parts of the 
refractive index, size (the refractive index is kept inde- 
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Fig. 2. Imaginary part of the refractive index for crys- 
talline forsterite (Mg 2 Si0 4 ) in the three different crystal- 
lographic directions. 



pendent of radiation wavelength), shape, structure. So far 
we have considered homogeneous and core-mantle spheres 
and infinite cylinders, and homogeneous spheroids, but are 
planning to involve more scatterer models. The figures in 
the DOP are given for weakly absorbing (n = 1.3 — 1.7, 
k = — 0.3) and highly absorbing (n = 1.5 — 3.5, 
k = 0.5 — 2.5) particles (see Fig. 3 for examples). These 
refractive indices are typical of astronomically interesting 
materials in visual. 

The special notes are the chapters of the review of 
Voshchinnikov (2002) which cover the following ques- 
tions of observed phenomena modelling — extinction 
efficiencies: general behaviour and deviations, wave- 
length dependence, the 2175 A feature, absolute ex- 
tinction and abundances, etc.; polarization efficiencies: 
size/shape/orientation effects, linear polarization: wave- 
length dependence, circular polarization: change of sign 
and so on. Some more chapters will be included in the 
near future. 
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Fig. 3. The extinction factors for core-mantle spheres (a) 
and infinite cylinders (b) of different core size (volume). 



3.5. Self-training algorithm for fractal-like aggregates 
(B5) 

This part of the DOP was developed in Kharkov 
University by V.P. Tishkovets, S.A. Beletsky and 
P.V. Litvinov. Its idea is as follows — as calculations of the 
optical properties of subparticle aggregates need a huge 
amount of processor time, it is worthwhile to construct 
a self-training algorithm that could learn to predict the 
properties analyzing a limited set of these data. 

To reduce the number of free parameters, fractal-like 
aggregates of equal size subparticles were selected. The 
structure of the aggregates can be characterized by two 
parameters: the fractal dimension D and the prefactor 
constant p (see, e.g., Feder, 1988). Other parameters are 
the number of subparticles N, the size parameter x and 
the real and imaginary parts of the refractive index of 
subparticles. 

The algorithm was based on an artificial neural net- 
work named 'multi-layered perceptron' (see, e.g., Miller, 
1990 for more details). A model of perceptron with 6 in- 
put and 30 output neurons and two layers with 30 neu- 
ron in each between the input and output ones was used. 
The algorithm was realized in the form of Pascal and 
Fortran codes. The region of its training was n — 1.4—1.7; 
k = 0.001 - 0.1; x = 1.5; D = 3; p = 8; N < 50. Accurate 
calculations were done by the Mackowski and Mishchenko 
(1996) code. The results were the expansion coefficients of 
the scattering matrix elements in series of the generalized 
spherical functions for particles of different N, refractive 
index, etc. Figure 4 shows the accuracy of the perceptron 
results. 

The codes and data archive as well as documentation 
are available via the DOP. Our work with the perceptron 
on modelling of the optical properties of interplanetary 
dust grains has demonstrated its advantages due to ex- 
tremely small size of the codes and very fast data access. 

4. Concluding remarks 

We have introduced a database of the optical properties 
of cosmic dust analogues. It is designed for student and 
scientists to help in applications of non-spherical scatterer 
models in astronomy but can be useful in other scientific 
fields as well. 

A further development of the database is planed in 
near future. It will include a graphics interface for the 
JPDOC, a special tool for on-line calculations of differ- 
ent optical properties of non-spherical layered particles, 
more chapters from Voshchinnikov's (2002) review and its 
continuation, and so on. 
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Fig. 4. Linear polarization degree P versus the scattering 
angle 9. Solid curve corresponds to the accurate values, 
dotted one to values given by the perceptron. The number 
of subparticles in an aggregate N — 20, m = 1.7 + 0.08i 
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